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Abstract 


An  analysis  is  presented  of  the  steady-state  propagation  of  a 
semi-infinite  mode  I  crack  for  an  infinite  inhomogeneous,  linearly 
viscoelastic  body.  The  shear  modulus  is  assumed  to  have  a  power-law 
dependence  on  depth  from  the  plane  of  the  crack.  Moreover,  both  a 
general  and  a  power-law  behavior  in  time  for  the  shear  modulus  are 
considered.  A  simple  closed  form  expression  for  the  normal  component 
of  stress  in  front  of  the  propagating  crack  is  derived  which  exhibits 
explicitly  the  form  of  the  stress  singularity  and  its  material  depend¬ 
ency.  The  crack  profile  is  examined  and  its  dependence  on  the  spatial 
and  time  behavior  of  the  shear  modulus  is  determined. 


;  Aceesslor; 

'  NTIS  g'a; 


1.  Introduction. 


In  this  paper  an  analysis  is  presented  of  the  quasi-static  problem 
of  a  semi-infinite  mode  I  crack  propagating  through  an  infinite,  inho¬ 
mogeneous,  isotropic  linearly  viscoelastic  body.  Viscoelastic  models 
involving  both  a  general  and  a  power-law  behavior  in  time  for  the 
shear  modulus  are  considered. 

One  of  the  few  investigations  of  a  boundary  value  problem  involv¬ 
ing  a  nonhomo geneo us  viscoelastic  body  was  carried  out  by  Walton  (1984) 
in  which  a  quasi-static  punch  problem  was  addressed  assuming  a  shear 
modulus  exhibiting  a  power-law  time  and  depth  dependence.  This  inves¬ 
tigation  was  carried  out  by  first  appealing  to  the  correspondence 
principle  and  then  reducing  the  boundary  value  problem  to  a  generalized 
Abel  integral  equation.  An  analysis  of  the  dynamic  steady  state  pro¬ 
pagation  of  a  semi-infinite  anti-plane  shear  crack  in  a  general  infi¬ 
nite  homogeneous  viscoelastic  body  was  presented  by  Walton  (1982) . 

Under  very  weak  assumptions  on  the  shear  modulus  the  problem  was  re¬ 
formulated  as  a  Riemann-Hilbert  problem  and  ultimately  simple  closed 
form  expressions  were  constructed  for  the  stress  intensity  factor  and 
the  entire  stress  field  in  front  of  the  advancing  crack.  The  approach 
adopted  in  this  study  employs  methods  from  both  papers.  The  corre¬ 
spondence  principle  is  utilized  to  derive  an  integral  relation  between 
the  normal  components  of  the  stress  and  displacement  fields  and  this 
relationship  is  subsequently  reduced  to  a  Riemann-Hilbert  problem. 

Before  proceeding  to  the  formulation  of  the  boundary  value  problem 
under  consideration,  a  few  comments  on  the  viscoelastic  models  assumed 
in  this  work  are  in  order.  The  homogeneous  power-law  viscoelastic 
model  may  be  characterized  by  a  shear  modulus,  y ,  of  the  form 


2 


u(t)  = 

where  t  denotes  time,  t  is  a  characteristic  relaxation  time,  and  y 

c  c 

is  a  characteristic  modulus.  Clearly,  this  simple  power-law  behavior 
in  time  introduces  certain  unphysical  features  into  the  model.  How¬ 
ever,  as  discussed  by  Walton  (1984)  and  Walton  et.al.  (1978)  amd  re¬ 
ferences  cited  therein,  the  pure  power-law  model  can  yield  useful 
engineering  predictions  for  many  real  materials. 

In  addition  to  the  power-law  model,  a  general  form  for  the  time 
dependence  of  the  shear  modulus  is  also  considered.  More  specifically, 
the  only  restrictions  placed  on  the  time  dependence  of  the  shear  mod¬ 
ulus  are  that  it  be  a  positive,  continuous,  decreasing,  and  convex 
function  of  time.  Though  the  solution  of  the  Riemann-Hilbert  problem 
requires  a  separate  analysis  depending  on  the  particular  form  of  the 
time  dependence  assumed  for  the  shear  modulus,  the  normal  component 
of  the  stress  field  in  front  of  the  crack  is  identical  in  both  cases. 
Indeed,  the  entire  stress  field  is  independent  of  any  time  dependence 
assumed  for  the  modulus. 

In  this  paper  the  viscoelastic  body  in  which  the  crack  propagates 
is  assumed  to  be  inhomogeneous  in  that  the  shear  modulus  has  a  power- 
law  depth  dependence.  In  particular  the  shear  modulus  y(t,y)  has  the 
form 

p(t,y)  =  y^,g{t)  (y/y^)^ 

where  y  denotes  distance  as  measured  from  the  plane  of  the  crack,  y^ 
is  a  characteristic  depth,  0^Y<1,  and  g(t)  satisfies  the  conditions 
discussed  above.  Again,  it  may  at  first  appear  that  a  model  with  a 
simple  power-law  depth  dependence  is  physically  unreasonable  and 
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therefore  only  of  mathematical  interest.  However,  as  pointed  out  by 
Walton  (1984)  there  are  materials  whose  mechanical  behavior  obeys  such 
a  power-law  model  over  very  broad  ranges  of  time  and  depth.  Though 
the  materials  and  examples  cited  concerning  a  power-law  depth  depend¬ 
ence  occur  more  naturally  within  the  realm  of  punch  problems,  there 
are  scenarios  for  which  the  power-law  model  applies  to  the  physical 
situations  that  arise  in  the  case  of  a  crack  in  a  viscoelastic  mater- 
i^l*  For  instance,  boundary  value  problems  involving  thermorheo- 
logically  simple  viscoelastic  materials  with  a  spatially  dependent 
but  time  independent  temperature  field  are  of  the  type  amenable  to 
analysis  under  isothermal  but  nonhomogeneous  spatial  assumptions. 
Within  this  context  the  power-law  model  may  provide  useful  insight 
when  the  mechanical  deformation  associated  with  the  propagation  of 
the  crack  produces  a  significant  thermal  gradient  which  in  turn  re¬ 
sults  in  a  depth-dependent  softening  of  the  material.  Other  physical 
situations  in  which  the  power-law  model  is  appropriate  occur  when  a 
crack  propagates  in  a  polymer  between  laminated  materials  where  mois¬ 
ture  has  lead  to  a  softening  of  the  polymer  material  and  in  certain 
instances  of  hydraulic  fracturing. 

In  the  next  section  the  transformation  of  the  governing  boundary 
value  problem  to  a  Riemann-Hilbert  problem  is  presented  and  the  solu¬ 
tion  of  the  Riemann-Hilbert  problem  is  derived.  In  Section  3  an  anal¬ 
ysis  of  the  stress  field  is  carried  out  and  it  is  shown  that  the  usual 
square-root  singularity  associated  with  the  stress  near  the  crack  tip 
is  altered  by  the  inhomogeneity  of  the  material.  A  very  simple  closed 
form  expression  for  the  normal  component  of  stress  is  derived  which 
exhibits  explicitly  the  nature  of  the  singularity  and  its  material 


dependency  and  indicates  that  the  usual  notion  of  the  stress  intensity 
factor  must  be  modified.  The  section  closes  with  an  asymptotic  analysis 
of  the  crack  profile  near  its  tip.  It  is  observed  that  the  displace¬ 
ments  of  the  crack  surfaces  depend  on  both  the  spatial  and  time  behavior 
of  the  shear  modulus. 

2.  Formulation  of  the  Problem. 


The  boundary  value  problem  to  be  studied  is  that  of  a  mode  I  semi 
infinite  crack  propagating  with  constant  speed  in  a  nonhomogeneous  iso 
tropic  viscoelastic  media  in  a  state  of  plane  strain.  Our  ultimate 
goal  is  to  describe  the  normal  component  of  stress  in  front  of  the 
advancing  crack  and  the  nature  of  the  crack  profile.  If  we  assume  the 
crack  propagates  in  a  plane  cdiout  which  the  spatial  properties  of  the 
body  are  symmetric  and  along  the  Xj^-axis  with  speed  v,  driven  by  loads 
f(Xj^-vt)  which  follow  it,  then  the  specific  boundary  value  problem  to 
be  solved  is 


j  =  0  -«=<Xj^<«>,X2>0 

a.  .  =  2u*de.  .  +  6.  .  (  2^)  y*de,  , 
1 :)  1  ]  1  ]  l-2v  kk 


O22ix^,0,t)  =  f(x^-vt)  x^  < vt 


U2(Xj^,0,t)  =0  Xj^  >  vt 


Oi2{Xi,0,t)  =  0  -  «>  <  x^  < 


Here  u^  are  the  viscoelastic  stress,  strain,  and  displace 

ment  fields,  v  is  Poisson's  ratio  (assumed  to  be  constant),  u  is  the 


shear  modulus,  6. .  is  the  Kronecker  delta,  and  u*de  denotes  the 


Riemann-Stielt jes  convolution 


vi*de 


•t 

u (t-T) de (t) . 

— oo 


The  construction  of  the  solution  of  the  preceding  equations  is  facili¬ 
tated  by  adoption  of  the  Galilecui  variable  x  =  Xj^-vt  and  the  change 

of  variables  =  ‘^xx'^12  “  ^xy''^!  ~  ’^x'  shear 

modulus,  two  cases  are  considered: 

case  (1)  Uj^(t,y)  =  u^g(t)  (y/y^)  ^ 

case  (2)  \i^{\L,y)  =  u^(t/t^)  “(y/y^)^. 

In  both  cases  0_<y  <1  and  in  case  (1)  g(t)  is  assumed  to  be  positive, 
continuously  differentiable,  nonincreasing,  convex,  and  such  that 
y  («>)  =  lim  y  (t)  >  0. 

A  key  step  towards  determining  and  u^  is  the  construction  of 
the  so-called  transfer  function,  i.e.,  the  function  T(p)  for  which 


°yy^P'0)  =  T(p)Uy(p,0) 


where  f(p,y)  denotes  the  Fourier  transform 


f(p,y)  = 


f (x,y)  e^^^dx. 


The  derivation  of  the  transfer  function  is  carried  out  as  in  Walton 
(1984) .  A  Fourier  transform  in  x  is  applied  to  the  boundary  value 
problem  and  then  we  appeal  to  the  viscoelastic  correspondence  prin¬ 
ciple.  Use  can  then  be  made  of  the  known  integral  boundary  relation 
between  elastic  stresses  and  displacements  for  a  power-law  depth- 
dependent  shear  modulus.  Specifically,  it  is  shown  in  Gladwell  (1980) 


(2) 


Uy(X.O) 


6  6  6 

where  u  ,  o  ,  and  a  are  the  elastic  displacements  and  stresses  euid 

y  yy 

sgn  ( • )  denotes  the  signum  function.  The  constants  g^2  922 

given  by 


g^2  =  (l--v)I  cos(7rq/2)/(iip®Y) 

922  =  (1-v)  Iq  sin  (•rrq/2)/i:Try®Y(l+Y) ) 
q  =  I(l+Y) (l-Yv/(l-v)) 

I  =  2^(Y+2)B(l±|±i,  1^). 


The  elastic  shear  modulus  adopted  by  Gladwell  has  the  form 

M®(y)  =  uV. 

The  viscoelastic  transfer  function  may  now  be  derived  by  Fourier  trans 

form  of  (2)  and  substitution  of  the  transformed  viscoelastic  shear  mod 
0 

ulus  for  (Recall  that  x  is  the  Galilean  variable  x,-vt.  Thus 

c  1 

ij(t,y)  is  a  function  of  x  and  y  and  m  denotes  its  Fourier  transform 
with  respect  to  x.)  If  use  is  made  of  (1),  it  can  be  shown  in  this 
way  that  the  transfer  functions  for  cases  (1)  and  (2)  are  as  follows: 

case  (1)  Oyy(p,0)  =  -K^ [ ( ipv) g (-vp) ] | p | ^“^Uy (p, 0)  (3) 

case  (2)  Syy(p,0)  =  -K^ [ ( vt^) “r (1-a) ( ip) j p | (p, 0)  (4) 

where 

U  r  (y+2)  cos{ytt/2) 

K  ^ 

o 


q(l-v)I  sin(qiT/2)y^ 


After  introduction  of  the  functions 


Hj^(p)  =  (ipv)g(-vp)  =  g{0)  + 


(5) 


and 


ia^  sgn(p) 

H2(p)  =  (vt^)  r(l-Qi)  {ip)“  =  (vt  )“r(l-a)  |pl“e 


(3)  and  (4)  may  be  written  as 


yy  yy 


(6) 


Here  and  denote  the  restriction  to  the  positive  and  negative 

real  axes  respectively  of  Similarly  for  u^  and  u^,.  Line  (6) 

may  be  viewed  as  the  Riemann-Hilbert  problem 


f'^(P)  =  -KqH^(p)  1p1^"Yf"(p)  +  h(p),  i=l,2 


(7) 


where 


F’^(P)  =  '*^,0)  =  ^ 


a  (x,0)e^P^dx, 

\  yy 


F  (p)  =  Uy(p,0)  =  ^ 


Uy (x,0) e^P^dx, 


and 


h(p)  =  -ayy(p,0)  = 


f (x) e^P^dx. 


It  is  assumed  a  priori  (and  it  may  be  verified  a  posteriori)  that 
Oyy  and  u^  are  such  that  and  F  define  functions  analytic  in  the 
upper  and  lower  half-planes  respectively.  Moreover,  the  limits 


F  (p)  =  lim  F  (p+iq),  F  (p)  =  lim  F  (p+iq) 
q-*-0  q->-0 


both  exist  and  satisfy  (7)  . 

Attention  is  now  directed  to  the  solution  of  the  Riemann-Hilbert 
problems  (7).  (The  reader  who  is  unfamiliar  with  the  theory  of  Riemann 
Hilbert  boundary  value  problems  may  wish  to  consult  Gaikov  (1966).)  To 


solve 


f'^(P)  =  -K^Hj^(p)  1p|^“^f"(p)  +  h(p) 


it  is  useful  to  consider  first  the  homogeneous  problem  of  finding 

+  ^  • 
functions  <I)  (z)  analytic  for  Imz<  0,  respectively,  which  satisfy  the 

homogeneous  boundary  relation 

^■^(p)  =  -K^Hj^(p)  1p|^"^<I>“(P)  .  (9) 

It  suffices  to  solve  separately 


(p)  =  (p)  '!>“  (p) 


=  -Kq  Ip  I 

—  +  + 

for  then  4>  (z)  =  ( z)  <I'2  ( z)  is  a  solution  to  (9).  In  order  to  solve 

(10)  it  is  necessary  to  determine  the  mapping  properties  of  Hj^(p). 
With  the  stated  assumptions  for  g(t) ,  it  is  shown  (Walton,  1982)  that 


i)  H^(0)  =  g(«>)  <  Re  Hj^(-vp)  £  g(0)  =  Hj^(“) 
ii)  ImH^(-vp)  =  -lmHj^(vp) 

(  >  0  p  >  0 


iii)  arg  H,  (-vp) 


From  these  properties  it  follows  that  the  image  of  in  the  complex 
plane  of  the  real  p-axis  is  as  depicted  in  Fig.  1.  In  particular 
Hj^(p)  is  Holder  continuous  for  all  p  on  the  extended  real  line  so 
that  the  canonical  solution  to  (10)  is  given  by 


(z)  =  exp  [j 


TTl 


T-Z 


dx] . 


Though  the  usual  methods  of  solving  a  Riemann-Hilbert  problem 
may  not  be  applied  to  (11) ,  since  the  coefficient  function  vanishes 
at  the  origin  and  becomes  infinite  at  <»,  a  solution  is  easily  found 
by  inspection.  Indeed,  if 


,+ ,  .  + 

$2(2)  =  c  z 


/1-Y. 


and  't'2(z)  =  c  z 


then 


<t2(P) 


c 

c 


(12) 


Hence  (12)  is  a  solution  to  (11)  provided  that  c^/c~  =  and 

^-{z)  =  c-z  expl^  J  ■  dT] 

—00 


is  therefore  a  solution  to  (9) . 

Before  proceeding  to  construct  a  solution  to  the  nonhomogeneous 
Riemann-Hilbert  problem  (8) ,  it  should  be  observed  that  ^^(z)  can  be 


greatly  simplified.  From  (5)  and  the  assumptions  on  g(t)  it  follows 
that  H, (p)  has  a  natural  extension  H, (p+iq)  for  q  <0.  Moreover 


i)  (±“+iq)  =  (“)  , 

ii)  <  Hj^(iq2)»  *32 


iii)  limH2^(iq)  =  Hj^(“) 
q ->•—<» 


Fig.  1  depicts  the  mapping  properties  of  Hj^(z)  for  Imz<  0.  Clearly 
for  Im  z  >  0 


,  f®  JlnH,  (t)  ,  ,  (•“  x,nM,  i  t; -x,nn. 

^  - - -  dT  =  i  JlnH,  (®)+^  - = - — — - 

TTl  T-2  2  1  2111  T-Z 


“  £,nH,  (  t)  (») 


From  the  stated  properties  of  H^(p-iq),  q  >0,  it  follows  that 
1  (  t) (®)  1^0  as  jx!  -*■<*>  cuid  hence  that  the  integral  in  (13) 

vanishes  by  an  application  of  Cauchy's  theorem.  Thus 

•f^Cz)  =  c^z  /g(  0) 


<I>‘’'(p)  =  c'*'/g(0)  k(p)  IpI 


where 


k(p)  = 


p  >  0 


iJd-Y) 


p  <  0 , 


The  solution  to  the  Riemann-Hilbert  problem  (8)  in  case  (1)  is 
now  determined  by  the  expression 


1  r  h(T)  dT 


Application  of  the  Plemelj  formulas  together  with  (14)  and  (15)  yields 

1-Y 


F*(p)  .  ^h(p)+  I 


h(T) 


dx 


1-Y  T-p  * 


(16) 


k(t)  It 


Attention  will  now  be  focused  upon  case  (2) ,  which  corresponds 
to  the  power-law  tinve  dependent  model.  To  solve 


f‘^(P)  =  -KqH2(p)  I  p|^"^F"(p)+h(p)  , 


(17) 


a  solution  to 


^■^(P)  =  -KqH2(p)  |p|^"^4'"(p) 


(18) 


+  ...+ 


is  sought  once  again  in  the  form  f  where 

’i'l(p)  »  -KQ(vt^)“rcl-a)  |p|^'^“"^'!'"(p) 


(19) 


cind 


'»'2(P)  =  e 


iay  sgn(p) 

•  Ilf 


»2(P). 

Problem  (19)  may  be  solved  as  before,  producing  the  solution 


(20) 


^^(z)  =  c*z 


where 


cVc'  =  -K  (vt  )“r(l-a) . 


(2 


A  solution  of  (20)  may  be  constructed  by  defining 


a(z)  =  S,n(z)  0  <  arg  z  <  2ir 


3(z)  =  iln(z)  -IT  <  arg  z<  ir, 
z  «  p+iq,  and  setting 

6*(z)  =  (a~(z)  +  e“(z))/2Tri  (22) 

+  + 

where  a  ,  6  denote  the  restrictions  of  a  and  8  to  the  upper  and  lower 
half-planes  respectively.  Then  6“(z)  is  analytic  for  Im  z <  0  respec¬ 
tively  and  it  can  be  verified  that 


6  (P) -6“(p)  =  -sgn(p) . 


Hence 


=  expt-ia  y  6*(z)l 


provides  a  solution  of  (20) .  Thus  a  solution  of  (18)  is  given  by 

±(Lt«7J.) 

H'‘(z)  =  c~z  ^  exp(-ia  j  6'(z)  1 


where  c^  emd  c  satisfy  (21)  and  6'*’  cuid  6~  are  defined  in  (22)  .  The 
calculation  of  ayy(p,0)  requires  that  ’J’'*^(p)  be  known,  euid  it  is  easy 
to  see  that 


l+g-Y  “ia  j  5^(p)  1 

'J''^(P)  =  c'^IpI  2  e 


iJd+a-Y) 


p  <  0 


l+a-1 


-a/2  1 


c'*’|p| 


p  >  0)  I  1 


p  <  0)  le 


ij(l+a-Y) 


p  >  0 


p  <  0 


=  c  k(p) IpI 


where  k(p)  is  defined  in  (15). 


Thus  the  solution  of  the  Riemamn-Hilbert  problem  (17)  is 


P*,Z)  =  1  f  ^  ^ 

2^1  J 


Again,  m2Jcing  use  of  the  Plemelj  formulas  and  (23)  it  follows  that 
for  case  (2),  F'*‘(p)  is  given  by  (16).  Since  F‘‘‘(p)  =  ^^(PrO)  »  it  is 
clear  that  for  both  cases  (1)  and  (2) , 

1-Y 


5yy(P.O)  =  ^h(p)  + 


r  h(T)  dT 

J  T-p* 


k(T) |t| 


(24) 


Denote  by  f''  the  inverse  Fourier  tremsform,  that  is 


f''(x) 


=  ^  I  e‘^'^f(p)dp. 


Since  h(p)  =  -0”y(p,O)  ,  Fourier  inversion  of  (24)  yields  for  x>0, 

1“Y 


^  1  \  — ^ 

mCO 

k(T)  |t| 


Again,  it  should  be  observed  that  (25)  gives  the  normal  component  of 
stress  in  front  of  the  crack  for  both  cases  (1)  and  (2) .  In  particu¬ 
lar  it  is  clear  that  the  stress  depends  only  on  the  spatial  properties 
exhibited  by  the  viscoelastic  shear  modulus.  In  the  next  section, 

(25)  is  written  in  a  very  single  form  which  explicitly  displays  the 
nature  of  the  singularity  of  near  the  crack  tip.  Futhermore, 

the  crack  profile  is  examined. 


Stress  and  Displacement  Analysis. 

In  order  to  cuialyze  x  >0  the  following  well  known 

identity  relating  the  Fourier  and  Hilbert  tr2msforms  is  recalled: 

H[f]  =  [-i  sgn  (t)  f''(T)  1 ''  (26) 


where 


Define  functions  a  and  6  by 


k(T)  sqn(T) 
1+Y 
i.iTT- 


B(T)  = 


k(T)  |t  I 


with  k(T)  given  by  (15) .  Prom  (26)  -  (28)  and  the  definitions  of 
h(p)  ,  the  expression  for  '^yy^*)  given  by  (25)  may  be  written  as 


,+  /..»  _  1  f*--ixp  h(p)  |j 

'yy'"*’  "  I?  J  ®  “^r 


1  d  f  o-i 

fF  d5?  J  ® 


-ixp  h(p) 


HI-Oyy31(p)  dp 


■li  sgn  (t)  (Oy,y*B)  (t)  ]  (p)  dp 


j  ^[S(p)  [sgn(T)  (o"  *B)  (t)  ]'tp)]''(x) 


i  d 


7  ^ta(T)  *  [sgn(T)  (a  *6)  (x)  ]  ]  (x) 


■  Ik  ”  n. 


.  •  • 


(Here  *  denotes  the  convolution  (f*g) (x)  *  J  f (x-x) g(T) dx.)  Prom 
(27)  a(x)  is  calculated  as 


a(x)  = 


^  f  >111)532^  e-i«dT 
|t| 


T=y 


2Tr  Ixl 


-if  (l-T) 

e  ^  (1+e  ^’'^1 


X  >0 
X  <  0 


Similarly  for  6(x), 


6(x) 


-ixt 


T=y 


dx 


k(x)  |t| 


X  >  0 
X  <  0 


(30) 


(31) 


By  utilizing  (30)  and  (31)  ,  interchcUiging  the  order  of  integration  in 
(29),  and  simplifying  the  coefficient,  the  stress  for  x  >0  may  be 
written  as 


(X) 


cos(^) 

2i 


(s) 


■  ■  dx 

I±y 


(x-x)  (x-s) 


ds. 


After  setting  x  -  s-(-t(x-s)  and  performing  the  differentiation,  the 
normal  stress  ahead  of  the  crack  is  easily  seen  to  be  given  by 


X  >  0 


(32) 


(X) 


cos(-i2-) 


TTX 


O^yTT? 


—00 


ds 

s-x' 


In  view  of  the  simple  form  of  (32)  certain  observations  are 

apparent.  By  setting  y=0  the  expression  for  the  known  stress  in 

front  of  a  uniformly  moving  crack  in  a  homogeneous  elastic  body  is 

+ 


recovered  (Willis,  1967).  Near  the  crack  tip,  as  x -►0, 


~^fi^YT72 


where 


K(y) 


cos 


tr 


lO  a  (s) 

ioC  - 


S 


ds. 


(33) 


The  effect  of  the  par6uneter  y  is  to  reduce  the  order  of  the  usual 
square  root  singularity  near  the  crack  tip.  This  seems  physically 
plausible  in  view  of  the  implications  of  a  Vcuiishing  shear  modulus 
in  the  crack  pl^me.  The  coefficient  K(y)  in  (33)  may  still  be  called 
a  stress  intensity  factor  in  the  sense  that  it  is  the  coefficient  of 
the  leading  singular  term  of  the  asymptotic  expemsion  of  the  stress 
at  the  crack  tip. 

To  illustrate  the  behavior  of  K(Y) »  consider  the  case  in  which 
the  crack  is  driven  by  constant  tractions  on  a  finite  interval  travel 
ing  with  the  crack  tip,  i.e. ,o”y (x,0) =-P [H (-x) -H (-a-x) ]  where  a  and  P 
are  constants  and  H  is  the  Heaviside  function.  Then 


We  conclude  this  section  by  examining  the  normal  crack  displace¬ 
ment.  The  calculation  of  the  displacement  of  the  crack  faces  is  ted¬ 
ious,  though  in  principle  it  is  easy  to  construct  an  expression  for 
Uy(x,0).  In  particular,  from  (6)  it  follows  that 

Uy(p,0)  =  K^^[H.(p)]"^|p|Y"^Syy(p,0) 

amd  consequently 

Uy(x,0)  =  k;^I(h“^(P)  (34) 


This  calculation  depends  heavily  on  the  explicit  form  of  the  time 
behavior  of  the  shear  modulus.  Even  for  the  power-law  model,  the 
Fourier  transform  in  (34)  can  only  be  explicitly  calculated  for  cer¬ 
tain  ranges  of  the  parcuneters  a  and  y,  A  more  efficient  method  of 
analyzing  the  crack  profile  and  deducing  the  effect  of  the  viscoelastic 
and  spatial  parameters  is  to  examine  the  derivative  of  the  normal 
displacement. 

Consider  first  the  power-law  behavior  in  time,  i.e.,  case  (2). 

From  (6)  and  (24)  it  follows  that 


Uy(p,0)  =  C^e 


-iaysgn(p)  (  -- 

yy 


l+a-Y 


iirlp 


f 

— 2 -  -o 


h(T) 


_ dx 

1-Y  T-p 


’k(T)  lx  I 


(35) 


where 


with  3  defined  in  (28). 

In  case  2(a),  ^  may  be  obtained  by  Fourier  inversion  of  (36). 

In  fact  for  this  case  it  is  possible  to  write 

Uj^(p)  = 

and 

U2(P)  =  s(p)  [i  sgn  (t)  (cFyy*6)  ] 

where 


r(x) 


X  >  0 


-C^r(l-a+Y) 

i 


1 

]  I 1-o+Y 
|x| 


/  T  »'  V 

cos(-ij-5 
cos  +a7r)  X  <  0 


and 


s(x) 


,3-2a+Y» 
^  ~'2  > 
3-2a+Y 


2Tr|xl  ^ 


0  X  >  0 

e^te^  +e  *  ]x<0 


Hence  one  obtains  the  following  expression  for  the  slope  to  the 
crack  displacement, 


u 


y»x 


(x,0) 


fO  _ 

j  Cyy  (u)  r(x-u)  du+ij  sgn(a))s(x-u))  B(a)-u)dti>du.  (38) 


9 

‘.v.v; 


It. 


To  discover  the  nature  of  the  crack  profile  at  its  tip  compute 

lim  u”  .  This  yields 
x-*>0 


u  (0 
y,x' 


)  = 


■C^r(l-a+Y)cos(^) 


V<“' 


(39) 


_ dh) 

-2a+Y) 
“2 - 


(w-u) 


T+ 

“7 


since 


du) 


r(n-Y-a)  r(.^-”;^~^) 


(3-2atY) 

°  U)  ^  (CJ-U) 


rFy 


V'“> 


u 


du 


(39)  may  be  simplified  to  show  u  (0  )  -  0.  Thus,  in  case  2(a),  the 

y  I X 

crack  surfaces  close  smoothly  to  form  a  cusp,  in  spite  of  the  presence 


2 


S 


of  a  singular  normal  stress. 

In  case  2(b)  the  crack  tip  is  no  longer  cusp-like  but  the  tan¬ 
gent  to  the  crack  tip  has  finite  slope.  It  is  apparent  that  for 
(l+2a-Y)/2=  1  (37)  is  still  valid  while 

-ioty  sgn(p)  ^ 

U2(p)  =0^6  k(p) sgn(p) H [-0yyBl (p) . 

In  this  case  Fourier  inversion  of  (36)  yields 


Uy (x,0) 


fO  _ 

Oyy (u) r(x-u) du  - 


(40) 


2ti' 


r  (a)  sin(a7r) 


n 


0) 


a 


duidudp . 


From  (40)  it  easily  follows  that  the  tangent  to  the  crack  tip  has 
finite  slope,  corresponding  to  the  crack  faces  meeting  at  a  sharp 
comer. 

The  case  of  a  homogeneous  viscoelastic  medium  with  a  power-law 
time  dependent  shear  modulus  corresponds  to  y  =  0  in  cases  2(a)  euid 
2(b).  In  case  2(a),  a  >1/2  and  the  crack  profile  is  cusp-like,  while 
for  2(b),  a  =  1/2  and  the  crack  surfaces  meet  at  a  sharp  corner. 

These  results  agree  with  the  findings  of  Walton  and  Nachman  (1979) . 

In  case  2(c)  the  temgent  to  the  crack  displacement  becomes  ver¬ 
tical  at  the  crack  tip.  Formulas  analogous  to  (38)  and  (39)  can  be 
derived  by  Fourier  inversion  of  (36).  However,  the  calculations  are 
lengthy.  If  an  asymptotic  analysis  of  u  is  carried  out,  the  nature 

y  ^  * 

of  the  crack  profile  is  more  easily  found. 


If  a  >Yf  e  L  (R)  ,  otherwise 


C  f(0)  -ia-  sgn(p) 

“l<P>  ~  |p|U<.-y'  ^ 


To  determine  the  asymptotic  behavior  of  U2  at  first  use  the  identity 

l/(T-p)  =  -l/p+T/(p(T-p) )  to  write 


H[-a"  e] 
yy 


'  iF  f-  *  «'P' 

““k(T)  I  T  I  ^ 


/3+y, 

where  G(p)  =  ^^(|p|  as  |p|-*-®-  From  this  it  follows  that 


U2(P)  ~  - 


-iai-  sgn  (p) 

C^e  ^  k(p) 


r 

L4-2a-Y 


J  iZX 

k(T)  ItI^ 


auid  hence 


U2(P) 


C  (•“> 
o 


h(T) 


k(T)  |t  1  I] 


.  TI 
-lotj 


iid+a-Y) 


p  — o 


From  (41)  ,  (42)  ,  and  the  Abeliein  theorems  for  the  Fourier  transform, 
it  follows  that  the  domineint  term  in  the  asymptotic  expansion  of 
for  X  near  0  is  given  by 


.  ir 

-ia>r 


„  i.«-(l+a-Y)  _ 


+^a-y 


+^a-y 


k  T  -[ 


The  purely  elastic  result  is  recovered  by  setting  a=0  in  (43) .  The 
homogeneous  case  is  then  found  by  letting  y  -*-0.  This  yields  that  for 
X  near  0, 


which  agrees  with  the  previously  known  result  (Willis,  1967) . 

From  (38),  (40),  and  (43)  it  is  apparent  that  the  crack  profile, 
unlike  the  stress  field,  depends  on  both  the  elastic  amd  viscoelastic 
parameters  y  and  a.  The  larger  the  value  of  y,  and  hence  the  softer 
the  material  near  the  crack  plane,  the  greater  the  crack  opening.  The 
effect  of  the  viscoelastic  parameter  a  is  to  reduce  the  crack  displace¬ 
ments.  In  particular,  when  the  viscoelastic  effects  are  dominant, 
i.e.,  cases  2(a)  and  2(b),  the  bltant  profile  associated  with  a  verti¬ 
cal  tangent  at  the  crack  tip  is  transformed  into  a  crack  opening  that 
is  cusp-like  or  has  finite  slope. 

One  final  remark  is  in  order.  The  simple  explicit  expressions 
for  u~  given  by  (38)  ,  (40),  emd  (43)  cannot  be  obtained  for  a  general 
time-dependent  shear  modulus.  However,  given  the  monotonicity  amd 
smoothness  assumptions  discussed  previously,  it  is  not  difficult  to 
adapt  the  methods  employed  for  the  pure  power-law  case  in  order  to 
prove  that  a  vertical  tangent  at  the  crack-tip  is  predicted  by  the 
more  general  models. 
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